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Abstract 

In the electroweak phase transition there arises the problem of baryon number 
washout by sphaleron transitions, which can be avoided if the phase transi- 
tion is strongly enough first order. The minimal supersymmetric standard 
model has just two Higgs doublets Hi and H 2 , while the next to minimal 
model, NMSSM, has an additional singlet, N, this latter giving rise to the 
helpful feature that the Higgs potential contains a tree level trilinear field 
term. We use the tunneling criterion for the existence of a first order elec- 
troweak phase change. A quantitative statistical analysis indicates that with 
parameters of the NMSSM satisfying the experimental constraints a strong 
first order phase change occurs in about 50% of cases. 
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1 Introduction 



It is well known that there is difficulty in sustaining the hypothesis JT|] of 
baryogenesis at the electroweak phase transition in the minimal standard 
model. In this standard model one Higgs doublet case the source of CP vi- 
olation is the CKM quark mixing matrix which is too small to explain the 
observed baryon to entropy ratio. Also, requirements on the phase transi- 
tion (in a perturbative treatment) seem to lead to a Higgs mass smaller than 
the experimental lower limit |^|, |3|. To overcome these difficulties, attention 
has been given to extensions of the minimal Standard Model, involving the 
addition of extra scalars [Q, |], ]5|, 0, [7L M. Prominent among these is the 
extension to two Higgs doublets; other rather natural such extensions, ne- 
cessitating two Higgs doublets, are supersymmetric theories. In the minimal 
super symmetric standard model, MSSM, the Higgs sector is just two dou- 
blets || [lO], [Hi], [12], |13| . As in the standard model the requirement on the 



electroweak phase transition is that it be first order so as to avoid washout of 
the baryon asymmetry immediately after the transition; it appears difficult 
to avoid this baryon washout in the MSSM [|T^ . 

Here we shall discuss the next-to-minimal model, NMSSM, which has 
additionally one singlet Higgs scalar [^|, ID]. Our treatment is perturbative; 



we do not treat the non-perturbative ideas which it has been suggested |14| 
might rehabilitate baryogenesis in the minimal, one Higgs doublet, Standard 
Model; and might when fully established provide a more secure basis for 
judgement on the various models. Additionally to the case of electroweak 
baryogenesis, it is important to examine the nature of the electroweak phase 
transition because of the many theories which do not create a net B-L in 
a phase change at higher energy; such theories are also menaced by baryon 
washout. 

In the absence of hard information we have to adopt a hypothesis on 
the SUSY breaking scale and on the spectrum of the particles. We follow 
a number of papers of recent years in taking the SUSY breaking scale, Ms, 
to be of the order of 1 TeV; we take perfect supersymmetry above that 
scale |n| [l(J. Then at Ms the quartic scalar couplings are fixed by the 
gauge couplings and two more parameters. We then use the renormalization 
group (RG) equations to run down the quartic couplings to the electroweak 
scale, where we investigate the nature of the phase change. There are also 
cubic and quadratic supersymmetry breaking couplings, and there results a 
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space of variable parameters in which we investigate what proportion leads 
to a strongly first order electroweak phase change, and so is compatible with 
electroweak baryogenesis. 

There has been quite considerable previous work on the electroweak phase 
change in the MSSM. We are not aware of so much on the NMSSM . The 
work of Pietroni |T7| has pointed up that the NMSSM, in contrast to the 
MSSM, has cubic terms in the scalar field potential at tree level leading to 
the possibility of a potential barrier in radial directions even at tree level; and 
that this diminishes the importance of the m 3 T term, in the high-T expansion 
of the T-dependent part, which has a vital role in most other theories of the 
phase change. That work uses a unitary gauge which we consider to be not 
so secure a basis for the consideration of phase changes as the Landau gauge 
which we use [§, [L8[ . 



A further difference is that we include the /i-term in the superpotential 
of the NMSSM; fi = is incorporated as a special case. Though there is 
a known naturalness problem with the magnitude of /i, the often preferred 
solution of setting \x equal to zero raises cosmic domain wall problems [|1^ 



When /i 7^ one often used |8], [L2| [L7]] definition of the critical temperature, 



being that temperature at which the curvature of the potential at the origin 
vanishes, can no longer be justified. 

We adopt the alternative definition, being that temperature when tun- 
neling first becomes possible from the high-T vacuum to the low-T vacuum 
having non-zero Higgs doublets expectation values. This criterion being well- 
known and discussed |T|, ^U), has on occasion actually been used in calcula- 



tions pi ], though it is more difficult to implement than the curvature crite- 



rion. This is an important change from most previous calculations, such as 



those of Ref. §, pf in the MSSM and Ref. |LJ in the NMSSM. Also this 
paper has the new feature that it compares the results from the two criteria 
when both can be implemented, that is in the fi = case. 



2 Formalism of the NMSSM 

The magnitudes of the parameters of the superpotential are significant for 
the electroweak phase transition. With the usual notation for the quark and 
lepton superfields, and summation over generations understood, the super- 
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potential is 



W = g u Qu c H 2 + g d Qd c H 1 + g e Le c H 



2 

+ 



+fiH 1 H 2 + \H X H 2 N - ^N 3 - rN (1) 

where Hi, H 2 are the Higgs doublet, and N the singlet, superfield 0. The cor- 
responding tree level scalar field potential, below the supersymmetry breaking 
scale, is 0, 0] 

V = ^(XxiHtHx) 2 + \ 2 {HlH 2 f) + 

(A 3 + \ a ){H{H x ){HIh 2 ) - A 4 \h\H 2 

(\ 5 HlHi + \ & H ] 2 H 2 )N*N + {\ 7 HiH 2 N* 2 + h.c.) + 
\ 8 (N*N) 2 + (\fi\ 2 + (\fj,*N + h.c.))(H\Hi + H\H 2 ) + 
m\HlHx + m 2 H ] 2 H 2 + mlN*N - (m A HiH 2 N + h.c.) - 

\(m 5 N 3 + h.c.) + -{m\H x H 2 + h.c.) + (m 2 N 2 + h.c.) (2) 

where 

H\ T = {Hi, Hi~), H 2 = (H 2 + , H 2 ), H X H 2 = Hi°H 2 ° — Hi H 2 + 

Hi, H 2 and A" now denoting pure scalar fields. The m-coefficient terms com- 
prise all possible soft supersymmetry breaking terms H |K]] . Vq is a function 
of 10 real scalar fields, the Higgs doublets and the singlet being given in 
terms of these by 

For simplicity, and to automatically ensure real VEVs, we shall follow the 
usual practice and take the parameters real. The boundary values at Ms of 
the quartic couplings are given by [p2fl 

A x = A 2 = |(sf + 9i), ^3 = ^(^2 - 9i), A 4 = A 2 - 



2 A possible term TV 2 in the superpotential can be removed by a field redefinition jf0| . 
The linear term, rN, contributes terms to the effective potential of the same type as we 
include with arbitrary coefficients in the soft supersymmetry-breaking part of the potential 
Vo given in Eq. (§). 
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Ak — As 



A 2 ,A 7 



-Xk, Ag 



and are developed down to M Weak by using the appropriate RG equations 
. The A and k are in principle free parameters at Ms @, TH ] . However they 
are linked to the one important Yukawa coupling [] , g t , by 3 simultaneous RG 
equations |23] in developing from high energy down to Ms; their values there 



should not be such that they correspond to divergent or unnaturally large 
values at high energy. From Eqs. (§) and (^|) we obtain the zero temperature 
potential Vo(0) as a function of (p = 02, ■ • ■ , 4>w with the parameters Aj 
assumed renormalised at the electroweak scale. The vacuum expectation 
values of the scalars are to be found at the minimum of Vq and thus satisfy 



diV n 



=<0> 



0. 



* 9 

Oi = 7T— , 



1,2, 



10 



(4) 



and also the requirement that (0) be a minimum. 

It is a constraint on our parameters that the minimum of Vq at zero 
temperature is for the VEVs of the Higgs fields having neutral components 
only 



TO 



AN) 



x 



(5) 



where V\, v 2 and x are real and v = Jv \ + v \ 
squared matrix is given by 



174 GeV. The scalar mass- 



M% = didjVo 



(6) 



and gives rise to 7 massive physical particles and 3 zero mass would-be Gold- 
stone bosons. The m\,m\^m\ are standard mass parameters and are to be 
specified in terms of the VEVs and the other parameters by Eq. (^) when (p 
is given by Eq. @. 

We can now discuss the other parameters in Vq. 

Firstly there are the terms involving /i which arise from the \x term in 
the superpotential. This raises the mu- problem (first noted in the MSSM) in 
the context of the NMSSM Jl9|, [2fJ. Unlike the A and k parameters of the 



superpotential just discussed there is no control on //, which would naturally 



3 we are not considering large tan (3 here 



4 



be expected to take on a value of the order of magnitude of the fundamental 
scale of the theory, whereas phenomenologically it should be of the order 
of the other electroweak terms. We do not take the point of view that the 
NMSSM can solve this by its having largely phenomenologically equivalent 
terms in the N field and simply setting \x = |23| , [25f . This can have its own 
difficulties when a resulting Z 3 symmetry gives rise to domain walls ||19|| . We 



tolerate the mu-problem. It should be noted that we have extra Z 3 breaking 
by the phenomenological term m\HiH2- Secondly there are the remaining 
soft parameters 777.4, ^5, 1^7', their provenance as completing the most general 



NMSSM breaking V was given in Refs. |g [TO] . 

Going now to the T-dependent terms in the effective potential we have 
the development, by the usual methods [H| for a temperature sufficiently 
high compared to the masses, given by 



V(<f>,T) = V (<f>) 



90 



+ V 2 (<f>,T) + V 3 (<f>,T) + V ln (<f>,T) (7) 



V 3 (<f>, T) 



T 



10 

i=i 



+ 4 



n, 



2 8 
^ i=l 



Mf_ 
64tt 2 



In 



+ 2 



Ml 

C;T 2 



i=l 



(8) 
(9) 



In Eq. (H) ± refers to fermions and bosons respectively, and is the number 
of degrees of freedom of the particle [El . 

V2 = -^Vtrace 



where V trace is the trace of the mass squared matrix 
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v t 



trace 



4m\ + 4m 2 + 2m i + 8 A* 2 + 



(3Ai + 2A 3 + A 4 + A 5 + -(3g 2 2 + ( ?1 2 ))(2/f 1 t J ff 1 ) + 

(3A 2 + 2A 3 + A 4 + A 6 + - A (?>gl + gf) + 3g 2 )(2HlH 2 ) + 
(2A 5 + 2A 6 + AX 8 )(2N*N) + 8 A// (AT + h.c.) + 



4|(/i + AiV)| 2 + 2A 2 (i7 1 



Ho 



') + 4k 2 N*N. 



(10) 



The last line is the 1-loop contribution of the Higgsinos in the case where 
the masses of the winos and binos are taken large, of order Ms, and is 
the sum of their squared masses. The size of the parameters gives these 
masses to be significantly less than Ms and thus they should be taken into 
account in the 1-loop corrections. In Eq. the first term is the contribution 
from the scalar masses. Here we have not just used the tree level masses 
squared, which are eigenvalues of didjVo as used in Eq. fllCf), but we use the 
eigenvalues of didj(Vo(</>) + 1^(0, T)). It has been shown f|, ^ that this is 
equivalent to correcting the one-loop potential by adding to the loop all ring 
diagrams ('daisy diagrams') which makes what would otherwise be imaginary 

3 

values of [M 2 ] 2 , for some values of </>, into real values JIB, [2(J. Also in Eq. 
(H), a correction to the gauge boson contribution has been included; in the 
contribution of the longitudinal polarisation excitations, there is suppression 
due to a temperature dependent 'Debye mass' factor. A simple treatment of 
this due to Dine et al. p is just to drop the longitudinal contribution, and 
we follow this prescription. 

A technical complication that arises is that if ^ > 2.2 for bosons or 
> 1.8 for fermions then the expansion of Eq. (Q) breaks down; this is 
important at T = T crit . We solve this by the method of Anderson and Hall 
U and replace the last three terms of Eq. (|?p by 



i (2vr) 2 V Mi 



15T 
1 H h 

8M 



(11) 



in the appropriate high 4^ region. Here n« is the number of degrees of freedom 
of particle i. 
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3 Phase Change Calculations 



Calculations with the kinetic equations for the dilution of baryonic charge 
just after the phase transition give a baryon preservation condition || |?J 

52=S1>{ (12) 

T C rit 

where £ varies between about 0.9 and 1.5 according to the gauge and other 
couplings of the theory. We take £ = 1. In practice in various theories just 
two criteria have been used to find the critical temperature, T cr i t , and authors 
have made a choice of either one or the other; we use both and compare. 
One is the temperature, T , at which, for decreasing T, the curvature of 
the effective potential V(<f), T) at (ft = (assumed to be the previous global 
minimum) first vanishes in the Higgs doublet neutral field directions. The 
other is the temperature, Tc, at which the value of V at a different minimum 
with a non-zero v' first becomes the global minimum of V(v[,v' 2 , x' , T); here 
the prime denotes general values of the field quantities as opposed to the 
specific values V\,v 2 ,x proper for the T = physical theory. As discussed 
below only for /x = is the first criterion, T crit = T , applicable. Now the 
shape of V(v [, v' 2 , x', T) depends on the theory parameters. It is the space 
of variable parameters which we search to assess in what part, and what 
proportion, of it the criterion Eq. (|l^) is satisfied. These variables are: 

1. [p,, X, k], occuring in the superpotential Eq. (HD; 

2. [tan/3 = v 2 /vi,x] (together with the known fixed value of v) specify 
the VEVs and replace m\, m\, m\; 

3. [Mch], the mass of the charged Higgs, is similarly given by an analytic 
formula, and we take it to replace as a variable; 

4. [m^,,mg,mj]; m\, when non-zero, breaks the Z3 symmetry even in the 
absence of \x. 

It is also necessary to consider the permissible range for the top Yukawa 
coupling, g t . We do not investigate the case where tan (3 is very large, this re- 
striction implying that gl 3> g\; there then are 3 simultaneous RG equations 
for X,k,g t P3| . These run from high energy down to M s , where the X,k,g t 
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should not be such that they correspond to divergent or unnaturally large 
values at high energy. The first limitation that this induces is gt(Ms) < 1.06 
which corresponds to gt{m t ) < 1.12. The experimental value of m t (m t ), be- 
ing greater than about 150 GeV, induces a lower limit on g t through the 
relation m t = g t vsm.(5, where v = 174 GeV, resulting in gt{fn t ) > 0.85. 

Within the above small available range for g t we select g t (M$) = 1 and 
we now discuss the relevant ranges of the above variables (i)-(iv): 

1. [p,, X,k]: We start by noting that the presence of /x, important fun- 
damentally, has an interesting consequence for the minimum of V at 
non-zero T. We discuss this in the high T expansion. As displayed 
above, Eq. (|10"D , Vi is a function of the fields with the only linear terms 
being those in (A/iiV + h.c). This means that for any non-zero T the 
origin 0; = cannot be a turning value and so cannot be a minimum 
for /i ^ 0. We are in a situation different from that in the MSSM but 



similar to that discussed by Choi and Volkas [21 in a different theory: 
the false vacuum is not at the origin, and the transition to the true 
vacuum is more complicated but amenable to investigation. It follows 
that only for the special case /i = is the critical temperature criterion 
of vanishing curvature at the origin, T crit = T , applicable. We have 
investigated the cases X/i = 0, ±20 GeV, ±40 GeV, ±60 GeV. The 
RG equations in A, k, g t also limit the ranges of A and k, which we find 
roughly require a/A 2 ± k 2 /g t to be less than about 0.8 f23|| . With our 
value g t (Ms) = 1 we have investigated the three cases (A, k) = (0.65, 
0.1); (0.1, 0.65); (0.5, 0.5). 

2. [tan j3, x] : The value of tan (3 is associated with that of the Yukawa cou- 
pling and the top mass through m t = g t v sin (3. We have investigated 
values of tan/3 in the range 1.6 to 3.0 in steps of 0.2 which, with our 
value of g t , corresponds to a running top mass of 150 to 173 GeV. The 
parameter x is rather free; the values considered are: x = O.lv, 0.5v, 
v, 2v. 

3. [Mg/J: Accepting the interpretation of the CLEO b — > 57 rate as re- 
quiring that M c h > 200 GeV, we have searched in the region 200 < 
M ch < 400 GeV. 
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4. [raj, 7715, mf] The range of searches of these variables has been 

my m 5 
" 6< (50G^F <6 ' - 6< (50G^F <6 ' _ ^ < JpOGeV) < 

In the ranges indicated above for M c h,m 5 ,ml and 777,7, we have conducted a 
coarse-grained search on a multidimensional grid over the parameter space. 

The first computing task is to restrict the parameter space to satisfy 
experimental constraints at T = 0. The conditions imposed are: 

(i) the masses are real, 

(ii) the minimum of V at the specified tan/3,?; (= 174 GeV) and X IS db 
global minimum in the space v[, v' 2 , x', 

(iii) the charginos and neutralinos have masses greater than 45 GeV and 30 
GeV respectively, and 

(iv) the mass of the lightest Higgs scalar is greater than 65 GeV. 

We call the set of these acceptable parameter sets the basis space. These 
conditions on the basis space are more restrictive than those imposed in our 
conference report |27|. This results in a much smaller basis space than in 



[27] and a much greater proportion of baryon preserving cases. 

The experimental limitation, (iv), on the lightest Higgs mass does not 
severely limit the basis space; it cuts out only about 1/3 of the sets which 
would otherwise pass into the basis space, whereas the first three conditions 
have a more severe effect. 



4 Results and Conclusions 

We need to find what proportion of the above basis space satisfies the baryon 
number preserving condition Eq. (JT2|). We shall take this as a measure of the 
ease of baryon preservation in the NMSSM theory. 

We have examined an initial grid of 105,840 parameter sets. Of these 
1,760, being 1.66%, pass through the tests (i)-(iv) of Section 3 for produc- 
ing an appropriate T = mass spectrum and form the basis space. We 
have then tested each parameter set in that space for whether it passes the 
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baryon preservation test, with the critical temperature defined as being that 
where tunneling from the higher temperature vacuum to the low temperature 
vacuum first becomes possible, as described in Section 3. 

The breakdown of these numbers into sectors of different /i, A, k and x, 
is given in Table 1. It is clear, both from the overall result and from the 
breakdown, that the condition of baryon number preservation just after the 
electroweak phase change does not impose a significant extra constraint on 
the parameters of the NMSSM from the present, statistical, point of view. 
Overall 50% of the basis space passes the test of baryon preservation. De- 
spite the large number of parameter sets in the initial grid, the sample of 
ultimately successful cases is around 900, too small to permit a meaningful 
delineation of the acceptable regions of our 9-dimensional parameter space. 
In the particular case of the parameters /i, A, k of the superpotential, Table 1 
gives an indication of the variation of the successful proportion with A and k; 
for /i the criteria select a basis space with a very strong bias towards negative 
(relative to x) and, after that selection, the baryon preservation condition 
gives a rather constant proportion over different values of fi and so does not 
select further. The 65 GeV lower bound we have placed on the mass of the 
lightest Higgs scalar is ultra-conservative, since it comes from analyses of 
experiments performed within the SM and the MSSM. In the NMSSM, with 
its extra mixing and different couplings, lower mass Higgs scalars could have 



escaped detection pq| . All parameter sets in our basis space yield a lightest 
Higgs scalar with a mass less than 120 GeV. The tendency is for lower mass 
lightest Higgs scalars to be associated with a greater proportion of baryon 
preserving cases (80% for a mass of 70 GeV but only 20% for 120 GeV). 

In Table 2 we show some results for the special case /i = 0. Of chief 
interest is the comparison between the zero curvature at ((f)) = criterion 
and the tunneling criterion for the critical temperature. Though there are 
some case by case differences the two criteria are statistically nearly the same; 
to our knowledge this is the first time such a detailed comparison has been 
made. It was found that the depth of the broken symmetry minimum relative 
to the value of V in the neighbourhood of the origin changed rapidly with 
temperature. As the perturbative effective potential formalism suffers from 



well-known problems near the origin particularly fll8| , Pq , [29|1 , the tunneling 
criterion is preferable to the frequently used curvature criterion, being less 
dependent on the value precisely at the origin. 

Our conclusions for the case [i = broadly agree with previous work on 
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this case JHJ, although our formalism differs, e. g. our gauge choice produces 
different T-dependent contributions to the scalar mass matrix. To treat the 
hitherto uninvestigated case of /i 7^ 0, we have had to change from the vanish- 
ing curvature criterion to the tunneling criterion for the critical temperature; 
we have found that values of /i up to about 100 GeV are acceptable, with 
well-chosen values of the other superpotential parameters (A, k). Unlike the 
Standard Model or the MSSM, the NMSSM can easily coexist with baryon 
number preservation. The essential difference is that cubic terms can already 
be present in the T = potential, so that the T-dependent terms modify but 
are not solely reponsible for the existence of a minimum away from the origin, 
v = 0, at the critical temperature. 
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Table 1: Showing the proportion of the basis space of parameter sets which 
satisfy the B( baryon)-preservation criterion, as a function of A, k and of the 
T = value of x (= (N)). In the basis space the T = criteria are satisfied. 
The values of /i covered are given by 60 GeV > A/i > —60 GeV. 



A k 


x/v 


B preserving 
%— age of basis 




0.1 


40 


0.65 0.1 


0.5 
1.0 


66 
82 




2.0 


59 


0.1 0.65 


1.0 


50 


0.5 0.5 


1.0 


26 



Table 2: /i=0: Comparing the results from the criterion of vanishing cur- 
vature of V at ((f)) = 0,T crit = T , with the tunneling criterion, T crit = T c . 
The baryon preserving proportion of the basis space and the average value 
of the critical temperature, T crit , are shown for both cases . 



A k x/v 


B preserving 
%— age for 

T C rit — Tq 


Average of 

Tcrit — T 


B preserving 
%— age for 

T C rit — Tc 


Average of 

Tcrit = T c 


0.65 0.1 1.0 


80 


70 GeV 


79 


71 GeV 


0.5 0.5 1.0 


16 


82 GeV 


11 


80 GeV 
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